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Abstract. - The integral relaxation time for the difference of numbers of particles in the two 
potential wells (IRT-A) for double-well classical spin systems is introduced. For the uniaxial 
symmetry, it is given by a quadrature. Unlike the previously introduced integral relaxation 
time for the magnetization, the IRT-A at low temperatures describes the rate of crossing the 
barrier under all conditions, including the strongly biased case. In the high-barrier case, the new 
integral relaxation time approaches the inverse of the lowest eigenvalue of the Fokker-Planck 
equation Ai. It can be more conveniently found by numerical methods than the latter. 



The rate of thermal activation Y of a dynamical system over a potential barrier AU 
between two wells is exponentially small in the low-temperature or high-barrier limit: T = 
Tq exp(— AU/T). In this limit, all definitions of the rate yield the same result. Analytical 
methods based on solving the Fokker-Planck equation (FPE) for the flow of the probability 
current over the barrier ware developed by Kramers Jl| for the particles (see ref. [|j for a 
review) and by Brown || for magnetic particles (classical spins) . In situations where analytical 
solution is problematic, the thermal activation rate T was usually identified with the lowest 
eigenvalue Ai of the Sturm-Liouville equation corresponding to the FPE and Ai was calculated 
numerically Q]. In the high-barrier case, Ai is exponentially small and describes the slow 
process of the exchange of particles between the two wells, whereas all other Aj, i = 2, 3, . .. 
are not exponentially small and correspond to intrawell processes. For spin problems in 
the high-temperature limit AU/T -C 1, the lowest eigenvalue Ai describes relaxation of the 
magnetization and it is, again, more important than other eigenvalues describing relaxation 
of higher moments of the distribution function. At intermediate temperatures AU /T ~ 1, the 
exclusive role of Ai in the overbarrier relaxation is no longer obvious, since other eigenvalues 
are of the same order as Ai and they can make contribution to this process, too. In addition, 
Ai cannot be calculated analytically in this temperature range even for the simplest models, 
such as that of uniaxially anisotropic spins with or without a longitudinal field. 

In fact, at intermediate temperatures the overbarrier escape rate is not unambiguously 
defined, and one can use various alternative definitions of T. One of the quantities measuring 
the escape rate is the mean first-passage time (MFPT) ||, H which for the simplest models can 
be calculated analytically in the whole range of temperatures. Another useful quantity which 
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is more appealing to physicists is the integral relaxation time (IRT) Tj nt which is defined by 
the area under the relaxation curve for the magnetization or polarization after an infinitesimal 
instantaneous change of the applied field |§| . Equivalently, the IRT is given by the coefficient 
in the linear term in the low-frequency expansion of the susceptibility. For the uniaxial 
spins in a longitudinal field, the low-frequency perturbation theory for the FPE leads to the 
analytical solution for the IRT in the whole temperature range. For more complicated models, 
such as a uniaxial spin in an oblique field, the low-frequency susceptibility can be obtained 
numerically from the matrix-continued- fraction solution of the FPE || . This routine is much 
more convenient than diagonalization of large matrices to obtain Ai . 

Numerical solutions for the r ln t and Ai performed in ref. [[)) for the uniaxial spin model in 
a longitudinal field revealed a striking feature. Whereas in zero field Tj~ t coincides with Ai 
in the whole temperature range to within only 1.2%, both quantities drastically deviate from 
each other at low temperatures for the fields exceeding some critical value. Whereas Ai is 
exponentially small, as it should be, tC/ in not exponentially small and increases with the 
increase of the strength of the potential relief. The physical explanation of this effect was given 
in ref. [ jTo| ]. In the biased low-temperature case, almost all particles or spins concentrate in the 
lower well, whereas the upper well is depleted. Then Tj nt captures the fast relaxation in the 
lower well (the intrawell process) rather than the slow overbarrier relaxation process. Since 
the intrawell process is not what one is usually interested in, the results of refs. |)|, |l0) put the 
usefulness of the concept of the integral relaxation time in a general case under question. 

The aim of this Letter is to introduce another type of the IRT which describes the process 
of the overbarrier relaxation in all situations. The key idea is to consider the responce of the 
population difference A = N + — N-, where N± are the numbers of particles in the two wells, 
to the changing bias field. Unlike the magnetization response used to calculate the IRT, the 
A response is insensitive to the redistribution of spins within the wells and always describes 
the overbarrier process. To emphasize the difference, it is convenient to designate the old 
(magnetization) IRT as IRT-M and the new (population) IRT as IRT- AT. 

Let us consider a single-domain magnetic particle with a classical Hamiltonian 

0H = v(n)-n-8£(t,(jj), 8£(t,w) = 5£ exp(-iwt), (1) 

where /3 = 1/T, n is the unit vector representing the direction of the particle's magnetization, 
5£(t,w) = pVM s SH(t,uj) is the reduced alternating field acting on the particle, V is the 
particle's volume and M s is the saturation magnetization. The all-temperature analytical 
solution for the integral relaxation time is possible if v(n) is axially symmetric and the 
alternating field is applied along the z axis. The mostly studied model is a uniaxial particle 
in the longitudinal field 

v(x) — — £x — era; 2 , x = cos9 (2) 

with 

£ ee (3VM S H, a ee PVM 2 S K, (3) 

where K > is the anisotropy constant creating a potential relief with two wells separated by 
a barrier. The top of the barrier corresponds to 

x = x rn = -h, h ee £/(2tr) = H/(2KM S ), (4) 

the barrier height Av = v(x m ) — v(—\) is given by 



Av = g - i + £ 2 /(4cx) = «t(1 - h) 2 , 



(5) 
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for h > 1 the barrier disappears. 

The dynamics of the classical single-domain magnetic particle is governed by the Fokker- 
Planck equation which was first formulated by Brown ||. This FPE can also be found in 
refs. H, [h| |ll). The equilibrium distribution function satisfying the static FPE can for axially 
symmetric models be written in the form 

= exp[-v(x)] ^ z= f 1 exp[ _ v(x)]dx _ (6) 



If an infinitesimally small alternating longitudinal field is applied, the distribution function 
can be searched in the form 

f(x,t,w) = f (x)[l + q(x)8t(t,wj\, (7) 
where the complex function q(x) satisfies the equation || [Io| 

— - — (1 - x 2 )^- + 2iuT N q = -(1 - x 2 )— - 2x. (8) 
ax ax J ax ax 

Here tm is the Neel relaxation time defined by 

=A N = 2 ia T/(VM s ), (9) 

where 7 is the gyromagnctic ratio and a is the dimcnsionless Gilbert's damping parameter. 
The solution of eq. ^ should satisfy the normalization condition J ^ fo(x)q(x)dx = 0. The 
low-frequency expansion of q{x) has the form 

q(x) = q (x) + iuiT N qi(x) + . . . , q Q (x) = x - B 

pX pi pX 

qi(x)= p(x')dx'- f (x) p(x')dx'dx, (10) 



where B — xfo(x)dx is the equilibrium spin polarization and 

p(x) = , 2 . $(x) = f (B - x')f Q (x')dx'. (11) 



l-x^f { X y 
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Having solved the FPE, one can calculate the dynamic susceptibility 

^ = lirwr = t x M*M*)te = x(i + + ...), x = B' = ^r. (12) 

Using eq. (|l(]) and integrating by parts one obtains the magnetization IRT [0 |[ |o) 

2t n f 1 dx $ 2 (x) 



B> J^l-x* f a (x) ■ 

In ref. @ it was shown that in the biased case £ > at low temperatures the integrand of this 
formula peaks near the top of the barrier x = x m and in the lower well x ~ 1 , which corresponds 
to the overbarrier and intrawell processes, respectively. For h > h c = 3 — 2v2 m 0.17 at low 
temperatures, the upper well is depleted, the overbarrier process is thus suppressed, and the 
intrawell process controls the IRT. This undesirable feature makes one to look for a more 
appropriate definition of the integral relaxation time. 
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Fig. 1. - Inverse of the integral relaxation time for populations IRT-iV, eq. and the lowest 
eigenvalue of the Fokker-Planck equation Ai || vs the potential strength a. 



The simplest choice for the quantity which is insensitive to the intrawell process is the 
difference of the wells populations 



A(i, w) = N + - iV_ = / sgn(x - x rn )f(x)dx. 



In equilibrium one has 

Z+ - Z- 



A = 



Z+ + Z_ 



2- = I fo(x)dx, 



Z+ = 



fo(x)dx, 



(14) 



(15) 



where Z± are the contributions to the partition function Z from both wells. The dynamic 
susceptibility for the population difference A has the form 



sgn(x - x m )f (x)q(x)dx = £a(1 + iuT iat ,N +■■■), (16) 



where 



- _ A/ _ 3A _ 2(Z_Z' + - Z + Z'_) 

XA ~ 9e (Z+ + Z.) 2 : 



and the integral relaxation time for the populations is given by 
2t n f 1 dx § N (x)$(x) 



Tint, W 



A' ./_! 1 - x 2 f (x) 



<f> N {x) = J [A - sgn(x' - x m )]f (x')dx' . 



(17) 



(18) 



This result follows from the same low- frequency solution of the FPE, eq. (10). It can be checked 
that the function &n(x) has no maxima near the bottom of the wells, thus Ti n t,JV does not 
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describe the intrawell process. This new integral relaxation time, which can be abbreviated 
as IRT- Y, captures the overbarrier process in the low-temperature range for all values of the 
longitudinal field. 

Dependences of Tj~ t N on the barrier-strength parameter a for different values of h for the 
uniaxial model [see eqs. (^) and (^)] are shown in Fig. 1, compared to the values of the lowest 
eigenvalue Ai quoted in Table I of ref . 0] . One can see that both quantities behave similarly, 
in contrast to r^ 1 (cf. Figs. 1 and 2 of ref. ||). In the region where the escape rate becomes 
exponentially small, r^ 1 N and Ai coincide with each other. An interesting feature of the 
dependences in Fig. 1 is the maximum of the relaxation rate arising in the strongly biased 
case, which is more pronounced for Ai than for r;~ t N . The physical reason for this maximum 
is the following. As long as the parameters a and £ are large but the reduced barrier height 
Av = <t(1 — h) 2 is small, the dynamics of the magnetic particle is qualitatively the same as 
that of a magnetic particle without a potential barrier at low temperatures. In this case the 
relaxation rate can be found in a simple way from the deterministic Landau-Lifshitz equation 
and it is proportional to the strength of the potential relief, i.e., to a. If er is further increased, 
then for any h < 1 the reduced barrier hight Av becomes large for a 3> 1/(1 — h) 2 , and the 
escape rate becomes exponentially small. These two regimes are separated by a maximum of 
the relaxation rate. 

In this Letter, the integral relaxation time for populations, IRT-Y, was introduced. This 
quantity is more useful than the previously considered integral relaxation time for the magne- 
tization, LRT-M, and it always captures the exponentially slow process of thermal activation 
in the high-barrier limit. For the simplest models considered here, there is an all-temperature 
solution for the IRT-Y in quadratures. In other cases, the problem can be solved numerically 
considering the low-frequency limit of the matrix-continued- fraction solution of the FPE . 
Formulation of the method for non-axially symmetric potentials or for (mechanical) particles 
poses no problems. Numerical calculation of the IRT-Y is a more convenient and accurate 
procedure than that for the lowest eigenvalue Ai. Using the new integral relaxation time should 
help to throw a light on fine details of various crossovers which have been recently discovered 
for the uniaxial model in a transverse field Jll| . For this model, numerical calculation of Ai 
faces difficulties in the low-damping high-barrier case. 

The author thanks Oleg Yevtushenko for critical reading of the manuscript. 
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